function ln_m = ln_m_KFE( m, ind, aux)
    % calculate grid for log(m) for distribution
    b = (log(m(2))-log(m(1)));
	ln_m = linspace(0,b,aux.n_m)';  
    if ind == 1
        ln_m = log(m(1)) + ln_m;
    end
    
end